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Abstract 

Matter-enhanced antineutrino flavor transformation between v e and v n( T ) can occur in 
supernovae if the vacuum masses for these species satisfy > m^ (r) . For Sm 2 > 1 eV 2 , 
such flavor transformation can affect the electron fraction Y e in the neutrino-heated super- 
nova ejecta. We point out that such flavor transformation will not drive Y e > 0.5 at the 
r-process nucleosynthesis epoch in the best available supernova model for such nucleosyn- 
thesis. Consequently, there is no obvious conflict between matter-enhanced antineutrino 
flavor transformation v e =f± v ^(t) an d r-process nucleosynthesis in the neutrino-heated 
supernova ejecta. 
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In this paper, we study matter-enhanced antineutrino flavor transformation between 
v e and v n( T ) m supernovae. In particular, we examine the effects of such flavor transfor- 
mation on the electron fraction Y e in the neutrino-heated supernova ejecta. We show that 
matter-enhanced flavor transformation v e =f± v ^(t) wn l n °t drive Y e > 0.5 when r-process 
nucleosynthesis takes place in the best available supernova model for such nucleosynthesis. 
In contrast to the case of matter-enhanced neutrino flavor transformation v e =f± v ^(t) ? there 
is no obvious conflict between antineutrino flavor transformation v e =f± v ^(t) an d r-process 
nucleosynthesis in the neutrino-heated supernova ejecta. 

The case of matter-enhanced neutrino flavor transformation v e =f± v ^(t) an d the resul- 
tant effects on r-process nucleosynthesis have been studied in Refs. [1] and [2]. The possible 
effects of matter-enhanced antineutrino flavor transformation v e =f± v ^(t) on r-process nu- 
cleosynthesis have been discussed in Ref. [3]. Reference [2] also discusses matter-enhanced 
antineutrino flavor transformation v e =f± v ^(t) m supernovae. It was pointed out in Ref. 
[1] that the electron fraction in the neutrino-heated supernova ejecta is determined by 
the characteristics of supernova neutrinos. In fact, the freeze-out value of Y e relevant for 
r-process nucleosynthesis is approximately given by the rates for the following reactions at 
the freeze-out radius: 

v e + n ^ p + e~, (la) 
v e + p =f± n + e + . (16) 

If we denote \ Ve n and \v eP as the rates for the forward reactions in Eqs. (la) and 
(lb), then the freeze-out value of Y e is given by 



In Eq. (2), the freeze-out radius rf Q satisfies 



K e n{rio) + Aj> eJ ,(r fo ) v(r {o )/r { , 



(3) 



where v(r{ Q ) is the outflow velocity of the ejecta at radius rf Q . The rates A I/en and Aj> eJ9 
decrease with increasing radius. Above the freeze-out radius, these rates become smaller 
than the rate at which the ejecta are flowing out. Consequently, v e and v e stop interacting 
effectively with the free nucleons in the ejecta above the freeze-out radius. 

In writing down the above equations, we have neglected the reverse reactions in Eqs. 
(la) and (lb). This is because the reverse reaction rates depend sensitively on the material 
temperature. At the freeze-out radius, the material temperature is low and these rates are 
small compared with the forward reaction rates. 

To calculate the rates \ Ve n and Aj> eJ9 , we need the angular and energy distributions 
of the supernova neutrino fluxes. Because the freeze-out radius is sufficiently large, we 
can make the approximation that supernova neutrinos are emitted from a neutrino sphere. 
This approximation gives a good description of the angular distribution of the neutrino 
fluxes at large radii. 

The cross sections for the forward reactions in Eqs. (la) and (lb) are given by 



respectively, where A np 1.293 MeV is the neutron-proton mass difference. The depen- 
dence of (T Ven on v e energy E Ve in Eq. (4a) is almost exact due to the Coulomb focusing 





and 




(46) 
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effect for the final state charged particles. The dependence of (Tjj eP on v e energy E„ e in 
Eq. (4b) is accurate when the positron in the final state is extremely relativistic. Because 
typical supernova neutrino energies are on the order of 10 MeV, this is almost always the 



case. 



From Eqs. (4a) and (4b), we note that analytic fits to the supernova neutrino energy 
distributions are appropriate for evaluating \ Ve n and \v eP if they can reproduce the first 
and second energy moments, (E„) and (E%), for individual neutrino flavors. In this case, 
the neutrino luminosity L v can serve as an overall normalization for the individual neutrino 
flux. The values of L„, (E„), and (E%) for v e , i> e , and are given in Table 1 for three 
different times relevant for supernova r-process nucleosynthesis: tpb (time post bounce) 
= 6, 10, and 16 s. The values of L„, (E„), and (E%) for v^, v T , and v T are approximately 
the same as those for v^. These values correspond to the neutrino characteristics obtained 
in the best available supernova model for r-process nucleosynthesis [4]. 

With these parameters, we can fit the local differential neutrino flux at radius r to 
the form: 

d<j> v ~ L \ \ A ; r dE v . (5) 

In Eq. (5), Fs(r] u ) is the rank 3 Fermi integral of argument r\ v . The values of the fitting 
parameters T v and T] v are also given in Table 1. 

From Table 1 and Eqs. (2), (4a), (4b), and (5), we can see that the freeze-out value 
of Y e is always less than 0.5 in the absence of neutrino flavor transformation. We give 
these values of Y e as in Table 1. The condition Y e < 0.5 is necessary for r-process 
nucleosynthesis to occur. We also note from Table 1 that the energy distributions for the 



v e , v e , and or neutrino fluxes are very different. As a result, the freeze-out value of Y e 
can be affected by transformation between different neutrino flavors. Since the v T and v T 
fluxes have identical energy distributions to the and fluxes, we will discuss the flavor 
transformation v e =f± and v e =f± as specific examples. It is to be understood that our 
conclusions regarding the effects of such flavor transformation on r-process nucleosynthesis 
are also valid for the flavor transformation v e =f± v t and v e =f± v t . 

References [1] and [2] show that matter-enhanced flavor transformation between a light 
v e and a with a cosmologically interesting mass (jn v 1 - 100 eV) can significantly 
change the freeze-out value of Y e . In fact, it is not necessary to have full conversion between 
v e and to drive Y e > 0.5 [1,2]. If r-process nucleosynthesis comes from the neutrino- 
heated supernova ejecta, then the required condition Y e < 0.5 places severe constraints on 
the mixing between v e and v^. 

However, the flavor transformation v e =f± can be enhanced in supernovae only if 
Vn is heavier than v e . So far there is no experimental evidence to confirm that this is 
the actual mass hierarchy. In fact, the current experimental upper limit on the v e mass 
is 7.2 eV [5]. Therefore, the possibility of v e and v e being the desired hot dark matter 
in some cosmological models [6] has not been ruled out yet. If the actual mass hierarchy 
is mjj e > m^, then the flavor transformation v e =f± can be enhanced in supernovae, 
whereas the flavor transformation v e =f± is suppressed. 

To illustrate the effects of the flavor transformation v e =f± i^, we plot the differential 
v e capture rate on protons d\/dE v = a „ eP dcf> v / ' dE ' v with respect to v e energy in Fig. 1. 
Since the differential capture rate is a radius-dependent quantity, the scale for the ordinate 
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is arbitrary. The solid line in Fig. 1 corresponds to the original v e energy distribution at 
tpb = 6 s, while the dashed line represents the case where the v e assume the same energy 
distribution as the at tpb = 6 s. 

From Fig. 1, we can see that the differential capture rate peaks at a neutrino energy of 
E v (E%) / (E„). This energy where the rate peaks is unique for each energy distribution. 
Between the peaks of the differential capture rates for the two energy distributions, the 
solid line and the dashed line cross at a neutrino energy of E v 25 MeV. Although Fig. 
1 is constructed for tpb = 6 s, we find that the generic features, especially the crossing of 
the differential capture rates for the two energy distributions at E v 25 MeV, as shown 
in Fig. 1 are also representative of the energy distributions throughout the period from 
tpb = 6 s to 16 s of the r-process nucleosynthesis epoch. 

It is obvious from Fig. 1 that the worst effect of the flavor transformation v e =f± 
on Y e occurs when antineutrinos with energies less than the crossing energy E v ~ 25 
MeV are fully transformed, whereas those with energies greater than the crossing energy 
remain unchanged. This is because there are more v e than with energies E v < 25 MeV. 
This worst case obtains, for example, when antineutrinos with energies E v < 25 MeV go 
through adiabatic flavor transformation and those with energies E v > 25 MeV do not have 
resonances below the freeze-out radius of Y e . This worst scenario may be realized because 
higher energy antineutrinos go through resonances at lower densities corresponding to 
larger radii. 

The resonance condition [7] for matter-enhanced antineutrino flavor transformation 

6 



=f± is given by 



Sm 2 



cos 26 = V2G F (n e + nf), 



(6a) 



where Sm 2 is the vacuum mass-squared difference, and 6 is the vacuum mixing angle. We 
choose Sm 2 > so that Sm 2 m? — m? for 6 <C 1. The net electron number density n e 
in Eq. (6a) is given by 



n e = n f 



n e + = Y e pN A , 



(66) 



where p is the matter density and Na is the Avogadro's number. 

The effective neutrino number density nf in Eq. (6a) represents the neutrino-neutrino 
forward scattering contributions to the neutrino propagation Hamiltonian [2]. It is defined 



to be 



n. 



n 



fAr)-nf(r)- nf(r)-nf(r) 



(6c) 



where, for example, nf(r) is the effective v e number density at radius r. 

In the case of matter-enhanced antineutrino flavor transformation, the flavor contents 
of the v e and fluxes are essentially unaffected for <C 1. With the approximation that 
neutrinos are emitted from a neutrino sphere at radius R v , the effective v e number density 
nf for a radially propagating neutrino is given by 



n. 



eff 



L, 



1 Rl 



(E Ve ) 47rr 2 c4r 2 



(6d) 



where c is the speed of light. The effective number density nf has an expression similar 
to Eq. (6d). The neutrino sphere radius at the r-process nucleosynthesis epoch is R v 10 
km. However, the effective v e and number densities nf and nf are affected by the 
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matter-enhanced antineutrino flavor transformation. They have to be determined from 
the flavor evolution history of individual antineutrinos with different energies. 

In the worst scenario given above, the freeze-out value of Y e will be larger than the 
original value Y" e a . However, for the neutrino energy distributions obtained in the best 
available model for r-process nucleosynthesis [4], we find that the freeze-out values of Y e 
in the worst scenario can never exceed 0.5 during the period from tpb = 6 s to 16 s when 
r-process nucleosynthesis takes place. We give these values of Y e as in Table 1. 

Of course, as discussed earlier, the worst scenario can happen only if antineutrinos 
with energies E v > 25 MeV do not go through resonances below the radius where the 
value of Y e freezes out. More specifically, the realization of this scenario depends on the 
matter density and effective neutrino number density at the freeze-out radius of Y e for a 
particular Sm 2 [cf. Eq. (6a)]. 

The freeze-out radius is approximately determined by Eq. (3). Since even in the worst 
scenario the total v e capture rate Aj> eJ9 is still greater than \ Ve n [cf. Eq. (2) and Table 1], 
we may conservatively approximate \ Ve n + ^v e p ~ 2A I/en . Together with Eqs. (4a) and (5), 
this gives 

ffoK r fo) ~ 10 14 cm 2 s -1 (7a) 

for the r-process nucleosynthesis epoch. 

One virtue of the best supernova model adopted for the r-process nucleosynthesis of 
Ref. [4] and in this paper is that it predicts the right amount of r-process material ejected 
from each supernova. In the context of galactic chemical evolution calculations the adopted 
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supernova model accounts nicely for the solar abundances of r-process elements [4]. The 
mass ejection rate M at radius r is given by 

M = 4:7rr 2 pv. (7b) 

This mass ejection rate stays at about 1O -6 M0 s _1 from tpb = 6 s to 16 s in the adopted 
supernova model [4] . 

We can rewrite Eqs. (7a) and (7b) as 

riop(rio) ~ 2 x 10 12 g cm -2 . (7c) 

From Eq. (7c), we can see that p(ri ) < 2 X 10 6 g cm -3 because the freeze-out radius 
always lies above the neutrino sphere, i.e., rf Q > R v ~ 10 km. In the actual supernova 
model, the freeze-out radius stays above 30 km. So the matter density at the freeze-out 
radius is below 7 X 10 5 g cm -3 . If we neglect the effective neutrino number density n e J^ 
in the resonance condition given by Eq. (6a), we find that antineutrinos with energies 
E v > 25 MeV will not go through resonances below the freeze-out radius for Sm 2 < 1 eV 2 . 

The effective neutrino number density is approximately 20 - 30% of the net electron 
number density at the freeze-out radius. As is evident from Eq. (6a), the effective neutrino 
number density tends to push the resonance position towards a larger radius for a given 
value of Sm 2 . This is to be contrasted with the opposite effects of the effective neutrino 
number density on matter-enhanced neutrino flavor transformation. As discussed in Ref. 
[2], for the case of matter-enhanced neutrino flavor transformation, the effective neutrino 
number density tends to draw the resonance position towards a smaller radius for a given 
Sm 2 . Taking into account the effective neutrino number density, we find that antineutrinos 
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with energies E v < 25 MeV will always have a resonance below the freeze-out radius if 
Sm 2 > 2 eV 2 . 

The scenario which gives the worst possible effect on Y e can then be realized for 
adiabatic antineutrino flavor transformation with Sm 2 2 eV 2 . For any other flavor 
evolution scenarios, the freeze-out value of Y e will always be smaller than in the worst 
scenario. It is interesting to observe that for a slightly larger 8m 2 3 eV 2 , antineutrinos 
with energies E v < 35 MeV will go through resonances below the freeze-out radius. If we 
further assume that these antineutrinos are fully transformed, then the freeze-out value of 
Y e is smaller than or very close to the original value Y" e a . We give these values of Y e as 
Y" e c in Table 1. For Sm 2 > 3 eV 2 , matter-enhanced antineutrino flavor transformation can 
decrease Y e below Y" e a . This is because at energies E v > 25 MeV, more are transformed 
into v e than v e are transformed into v^. This effect can be seen from Fig. 1. 

In the above discussion, we only have considered the possible effects of matter- 
enhanced antineutrino flavor transformation v e =f± on r-process nucleosynthesis in super- 
novae. However, if the mass hierarchy > obtains, there will be other consequences 
of the flavor transformation v e =f± v^. In particular, the flavor transformation v e =f± 
taking place in a galactic supernova could cause observable effects in future neutrino de- 
tectors such as super Kamiokande. This is because the flavor transformation v e =f± can 
affect the energy distribution of v e detected on earth. 

In fact, the predecessor of super Kamiokande, the Kamiokande II detector, detected 
11 neutrino events from SN1987A, while 8 events were detected by the now inactive 1MB 
detector. There have been many studies in the literature which try to extract information 
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about supernova neutrinos from these 19 (mainly v e ) events. Before we give a brief dis- 
cussion of these previous studies, we emphasize that the period of the supernova process 
relevant for supernova neutrino detection is much earlier than the r-process nucleosynthesis 
epoch. The majority of the neutrino events in a detector come from the first a few seconds 
after core bounce (i.e., tpb < 2 s), whereas the r-process nucleosynthesis takes place at 
tpb > 6 s. Both the neutrino characteristics (e.g., luminosities and energy distributions) 
and the dynamic aspects of the supernova (e.g., density structure and hydrodynamic in- 
stabilities) are very different for these two periods. The discrepancy between (E^ e ) and 
{E& ) at earlier times is smaller than as shown in Table 1 for tpb > 6 s. The neutrino 
sphere radius shrinks significantly from R v 50 km at tpb ~ 0.1 s to R v 10 km at tpb 
> 3 s. Consequently, treatment of the flavor transformation v e =f± and its implications 
will also be different. 

A careful statistical analysis of the SN1987A neutrino signals has been carried out in 
Ref. [8]. However, the results presented as Figs. 5 and 6 in Ref. [8] were based on an overly- 
simplified model of supernova neutrino emission. To be specific, it was assumed that (1) the 
v e energy distribution is given by a Fermi-Dirac distribution with zero chemical potential; 
(2) the temperature characterizing the Fermi-Dirac distribution decreases exponentially 
with time; (3) the neutrino sphere radius is fixed; and (4) the v e luminosity is given by 
the blackbody radiation law corrected for the Fermi-Dirac statistics, i.e., there exists a 
specific relation between neutrino luminosity, neutrino temperature, and neutrino sphere 
radius. These assumptions differ from the approximations for supernova neutrino fluxes 
made in this paper. While we employ a Fermi-Dirac distribution with finite chemical 
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potential for the normalized neutrino energy distribution, we treat the neutrino luminosity 
as an independent quantity [see Eq. (5)]. There is no simple analytic function which can 
adequately describe the evolution of neutrino luminosity and neutrino energy distributions 
simultaneously. In addtion, the above assumption (2) contradicts the physical effects of 
neutronization on the v e opacities. Protons in the core provide an important source for 
the v e opacity through the forward reaction in Eq. (lb). The lower opacity then hardens 
the v e energy distribution as the core becomes more deficient in protons with time. 

Given the poor statistics of the SN1987A neutrino events, and the above mentioned 
errors in the supernova neutrino emission model used to analyze these events, it is clear 
that conclusions regarding the v e energy distribution so obtained must be taken with great 
caution. With a total of only 19 events, one can hardly expect more than a confirmation 
of the gross energetics and qualitative features of supernova neutrino emission. With 
the above mentioned caveats on the underlying supernova neutrino emission model used 
to analyze the SN1987A data, fits to the v e temperature range between ~ 3 and ~ 6.5 
MeV [8]. This range is clearly too broad to preclude the possible occurrence of flavor 
transformation v e =f± in SN1987A. 

It is our sincere hope that nature will grant us the opportunity to detect neutrinos 
from a galactic supernova when super Kamiokande is in full operation. In that case, we 
will have many more neutrino events, which may enable us to extract details of supernova 
neutrino energy distributions. 

In conclusion, we have studied the effects of matter-enhanced antineutrino flavor 
transformation on the freeze-out value of Y e in the neutrino-heated supernova ejecta. We 
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find that such flavor transformation can never drive Y e > 0.5 when r-process nucleosyn- 
thesis takes place in the best available supernova model for such nucleosynthesis. For 
Sm 2 > 3 eV 2 , matter-enhanced antineutrino flavor transformation can even decrease Y e 
below the original value Y" e a for no flavor transformation. While the actual effects of dif- 
ferent Y e from Y" e a on the r-process nucleosynthesis still await further detailed study, there 
is no obvious conflict between matter-enhanced antineutrino flavor transformation and 
r-process nucleosynthesis in the neutrino-heated supernova ejecta. 

We acknowledge Jim Wilson for providing us with the details of his numerical su- 
pernova model used for r-process nucleosynthesis. This work was supported by the De- 
partment of Energy under Grant No. DE-FG06-90ER40561 at the Institute for Nuclear 
Theory and by NSF Grant No. PHY-9121623 and an IGPP minigrant at UCSD. 
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Figure Caption 

Fig. 1 Differential v e capture rate on proton with respect to v e energy. The solid line 
corresponds to the original v e energy distribution at tpb = 6 s. The dashed line obtains 
when v e assume the same energy distribution as at tpb = 6 s. 
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Table 1: Supernova neutrino characteristics and freeze-out values of Y e for different 
scenarios of antineutrino flavor transformation at three representative times during 
the r-process nucleosynthesis epoch. 



tpb (s) 


V 


L v 

(10 51 erg s" 1 ) 


(MeV) 


{El) 
(MeV 2 ) 


T v 
(MeV) 




Y a 

e 






6 


V e 


0.98 


10.5 


133.5 


2.69 


2.78 


0.44 


0.48 


0.42 




1.05 


18.4 


389.6 


3.69 


4.97 


v» 


1.64 


26.1 


901.1 


8.56 


-1.27 


10 




0.57 


10.2 


125.1 


2.55 


3.03 


0.40 


0.41 


0.34 




0.64 


19.5 


440.6 


4.08 


4.60 


v» 


1.26 


26.2 


915.0 


8.73 


-4.72 
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0.49 


10.1 


124.4 


2.66 


2.53 


0.37 


0.41 


0.37 




0.62 


19.6 


449.2 


4.31 


4.16 




0.90 


26.6 


937.4 


8.75 


-1.51 
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